home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
The Datafile PD-CD 1 Issue 2
/
PDCD-1 - Issue 02.iso
/
_utilities
/
utilities
/
001
/
meschach
/
!Meschach
/
c
/
matlab
< prev
next >
Wrap
Text File
|
1994-01-13
|
5KB
|
197 lines
/**************************************************************************
**
** Copyright (C) 1993 David E. Steward & Zbigniew Leyk, all rights reserved.
**
** Meschach Library
**
** This Meschach Library is provided "as is" without any express
** or implied warranty of any kind with respect to this software.
** In particular the authors shall not be liable for any direct,
** indirect, special, incidental or consequential damages arising
** in any way from use of the software.
**
** Everyone is granted permission to copy, modify and redistribute this
** Meschach Library, provided:
** 1. All copies contain this copyright notice.
** 2. All modified copies shall carry a notice stating who
** made the last modification and the date of such modification.
** 3. No charge is made for this software or works derived from it.
** This clause shall not be construed as constraining other software
** distributed on the same medium as this software, nor is a
** distribution fee considered a charge.
**
***************************************************************************/
/*
This file contains routines for import/exporting data to/from
MATLAB. The main routines are:
MAT *m_save(FILE *fp,MAT *A,char *name)
VEC *v_save(FILE *fp,VEC *x,char *name)
MAT *m_load(FILE *fp,char **name)
*/
#include <stdio.h>
#include "matrix.h"
#include "matlab.h"
static char rcsid[] = "$Id: matlab.c,v 1.7 1994/01/13 05:30:09 des Exp $";
/* m_save -- save matrix in ".mat" file for MATLAB
-- returns matrix to be saved */
MAT *m_save(fp,A,name)
FILE *fp;
MAT *A;
char *name;
{
int i;
matlab mat;
if ( ! A )
error(E_NULL,"m_save");
mat.type = 1000*MACH_ID + 100*ORDER + 10*PRECISION + 0;
mat.m = A->m;
mat.n = A->n;
mat.imag = FALSE;
mat.namlen = (name == (char *)NULL) ? 1 : strlen(name)+1;
/* write header */
fwrite(&mat,sizeof(matlab),1,fp);
/* write name */
if ( name == (char *)NULL )
fwrite("",sizeof(char),1,fp);
else
fwrite(name,sizeof(char),(int)(mat.namlen),fp);
/* write actual data */
for ( i = 0; i < A->m; i++ )
fwrite(A->me[i],sizeof(Real),(int)(A->n),fp);
return A;
}
/* v_save -- save vector in ".mat" file for MATLAB
-- saves it as a row vector
-- returns vector to be saved */
VEC *v_save(fp,x,name)
FILE *fp;
VEC *x;
char *name;
{
matlab mat;
if ( ! x )
error(E_NULL,"v_save");
mat.type = 1000*MACH_ID + 100*ORDER + 10*PRECISION + 0;
mat.m = x->dim;
mat.n = 1;
mat.imag = FALSE;
mat.namlen = (name == (char *)NULL) ? 1 : strlen(name)+1;
/* write header */
fwrite(&mat,sizeof(matlab),1,fp);
/* write name */
if ( name == (char *)NULL )
fwrite("",sizeof(char),1,fp);
else
fwrite(name,sizeof(char),(int)(mat.namlen),fp);
/* write actual data */
fwrite(x->ve,sizeof(Real),(int)(x->dim),fp);
return x;
}
/* d_save -- save double in ".mat" file for MATLAB
-- saves it as a row vector
-- returns vector to be saved */
double d_save(fp,x,name)
FILE *fp;
double x;
char *name;
{
matlab mat;
Real x1 = x;
mat.type = 1000*MACH_ID + 100*ORDER + 10*PRECISION + 0;
mat.m = 1;
mat.n = 1;
mat.imag = FALSE;
mat.namlen = (name == (char *)NULL) ? 1 : strlen(name)+1;
/* write header */
fwrite(&mat,sizeof(matlab),1,fp);
/* write name */
if ( name == (char *)NULL )
fwrite("",sizeof(char),1,fp);
else
fwrite(name,sizeof(char),(int)(mat.namlen),fp);
/* write actual data */
fwrite(&x1,sizeof(Real),1,fp);
return x;
}
/* m_load -- loads in a ".mat" file variable as produced by MATLAB
-- matrix returned; imaginary parts ignored */
MAT *m_load(fp,name)
FILE *fp;
char **name;
{
MAT *A;
int i;
int m_flag, o_flag, p_flag, t_flag;
float f_temp;
Real d_temp;
matlab mat;
if ( fread(&mat,sizeof(matlab),1,fp) != 1 )
error(E_FORMAT,"m_load");
if ( mat.type >= 10000 ) /* don't load a sparse matrix! */
error(E_FORMAT,"m_load");
m_flag = (mat.type/1000) % 10;
o_flag = (mat.type/100) % 10;
p_flag = (mat.type/10) % 10;
t_flag = (mat.type) % 10;
if ( m_flag != MACH_ID )
error(E_FORMAT,"m_load");
if ( t_flag != 0 )
error(E_FORMAT,"m_load");
if ( p_flag != DOUBLE_PREC && p_flag != SINGLE_PREC )
error(E_FORMAT,"m_load");
*name = (char *)malloc((unsigned)(mat.namlen)+1);
if ( fread(*name,sizeof(char),(unsigned)(mat.namlen),fp) == 0 )
error(E_FORMAT,"m_load");
A = m_get((unsigned)(mat.m),(unsigned)(mat.n));
for ( i = 0; i < A->m*A->n; i++ )
{
if ( p_flag == DOUBLE_PREC )
fread(&d_temp,sizeof(double),1,fp);
else
{
fread(&f_temp,sizeof(float),1,fp);
d_temp = f_temp;
}
if ( o_flag == ROW_ORDER )
A->me[i / A->n][i % A->n] = d_temp;
else if ( o_flag == COL_ORDER )
A->me[i % A->m][i / A->m] = d_temp;
else
error(E_FORMAT,"m_load");
}
if ( mat.imag ) /* skip imaginary part */
for ( i = 0; i < A->m*A->n; i++ )
{
if ( p_flag == DOUBLE_PREC )
fread(&d_temp,sizeof(double),1,fp);
else
fread(&f_temp,sizeof(float),1,fp);
}
return A;
}